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Abstract 

Taking the reported high-Q 2 anomaly at HERA as a signal of new physics, 
we show that several independent considerations point towards s-channel lep- 
toquark production as the most attractive interpretation. We argue that even 
this option is highly constrained by flavor-changing processes: the couplings 
must be accurately diagonal in the quark and lepton mass eigenstate basis and 
should preserve individual quark and lepton family numbers. We propose a 
dynamical mechanism that might produce this pattern; it has distinctive ex- 
perimental consequences. 
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Recently, the HI and ZEUS || collaborations, which have been studying e ± p 
collisions on the HERA ring since 1993, announced an anomaly in high-Q 2 e + p colli- 
sions. Using a combined accumulated luminosity of 34.3 pb -1 in e + p — > eX mode at 
y/s = 300 GeV, the two experiments have observed 24 events with Q 2 > 15000 GeV 2 
against a Standard Model (SM) expectation of 13.4 ± 1.0, and 6 events with Q 2 > 
25000 GeV 2 against an expectation of only 1.52±0.18. If we take these reports at face 
value, there seems to be little possibility of explaining the anomaly within the Stan- 
dard Model. In particular, the high-Q 2 events are clustered at Bjorken-x values near 
0.4 to 0.5; the quark distribution functions at such large x are well-measured at lower 
Q 2 , and QCD allows for a very precise scaling to the observed Q 2 . (By contrast, the 
superficially similar reported anomaly in the high-.E/r dijet cross-section at FNAL || 
can plausibly be ascribed to uncertainties in the gluon distribution functions.) The 
peaked distribution of these high-Q 2 events in x is also notable. If the new physics 
is assumed to be an s-channel resonance, this distribution translates into a mass de- 
termination. For both HI and ZEUS two independent determinations are possible, 
depending on whether x is calculated from the double-angle or electron methods. 
For the seven events selected for special study by HI we find these two methods give 
M 2a = 202 ± 14 GeV and M e = 199 ± 8 GeV respectively. Similarly, for the five events 
selected by ZEUS, M 2a = 231 ± 16 GeV and M e = 219 ± 12 GeV . 

HERA is capable of running in two modes: e~p and e + p. In the former mode 
HI and ZEUS have accumulated 1.53 pb -1 of data but have observed no statistically 
significant deviations from the SM. Further, the experiments differentiate between 
final state eX and uX, where the neutrino is detected through its missing p?. The 
anomaly discussed above is in the neutral current (NC) reaction e + p — > eX channel. 
HI has also announced its findings in the e + p — *> vX charged current (CC) channel. 
They find 3 events at Q 2 > 20000 GeV 2 with an expectation of 0.74 ± 0.39, but no 
events with Q 2 > 25000 GeV 2 . ZEUS has not announced its CC data. Compared to 
the NC channel, the CC signal is considerably less statistically significant, and we 
therefore regard the question of whether new physics in this channel is required as 
open. 

Both detectors at HERA are highly efficient at tagging muons. Apart from one 
such event in the 1994 HI data [|J (which occurs in a very different kinematic region 
from the events under consideration here), neither experiment has evidence for final 
state muons. 

Reprising: HI and ZEUS see an excess of high-Q 2 NC events in 34.3 pb _1 of e + p 
data, with no corresponding anomaly in 1.53 pb^ 1 of e~p data. There is no clear 
evidence for or against a similar anomaly in the CC mode. There is no signal for 
final state /z's. There is evidence that the events seen are peaked at x ~ 0.45-0.55, 
but there is some disagreement between the two experiments concerning the central 
x value at the peak. If the signal is real, and not a statistical fluctuation, there 
appears to be no way to change the parameters of the SM to explain it. Let us 
entertain the hypothesis that this data represents a genuine signal, and investigate 
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the consequences. 

New Physics. What kind of new physics could be responsible for the HERA 
anomaly? We consider only physics that can be described by an SU(3) x SU(2) x £7(1)- 
invariant effective Lagrangian using perturbation theory^ 

There are 3 kinematic channels through which the ep scattering might be occur- 
ring: in the s-channel through a resonance with non-zero lepton and baryon numbers 
(a "leptoquark"), in the w-channel again via a leptoquark, or in the t-channel via 
a particle without lepton or baryon number (this includes Higgs or Z' exchange or 
contact interactions). We will discuss the t and w-channels first. 

If the intermediate state is a t-channel scalar or vector, then there are two kine- 
matic limits in which very different constraints play a role. For a light Z' (the results 
hold equally well for a scalar with a suitably rescaled coupling), the CDF and DO 
experiments provide strong bounds from searches for qq — > Z' — ► e + e~,fi + fi~. As- 
suming significant branching ratios of the Z' to leptons, they rule out Z' masses below 
about 650 GeV |§ . For masses greater than this bound, the Z' can be integrated out 
and one can work in terms of effective 4-fermion contact interactions. By convention, 
one writes the operators: 

C = 1T¥T2 rfij {ea^i ) (QjYQj ) ( i ) 

i,j=L,R VHj) 

where q = u,d and rj^ = ±1. Thus there are 8 different contact interactions, times 2 
signs of rj. 

To analyze how well each operator fits the observed data at HERA, we perform 
a standard \ 2 analysis with one free parameter, namely Afj. For the fit, we use the 
combined bin-by-bin data of ZEUS and HI for bins starting at Q 2 =(10000, 15000, 
20000, 25000, 30000, 35000) GeV 2 respectively, with each bin (except the last) being 
5000 GeV 2 wide; the last bin goes from Q 2 = 35000 GeV 2 to Q 2 = s = 90000 GeV 2 . 
Bin by bin, ZEUS and HI combined see (29, 14, 4, 2, 2, and 2) events. 

For each operator and each choice of Af- we calculated the number of events 
expected in each Q 2 bin at HERA. HI and ZEUS claim efficiences approaching 80% 
in the Q 2 region of interest, and our own calculations of the SM and new physics cross- 
sections were scaled appropriately. For the parton-level calculation, we generalized 
the results of Ref. [|7| to contact interactions. The parton-level calculations were then 
folded into the CTEQ3M parton distribution functions || to get the physical event 
rates. 

The Standard Model, with 6 degrees of freedom has a x 2 /^o/=4.0, which is dis- 
favored at greater than 99.9% C.L. For the contact terms, we demand a fit better 
than 2cr (95% C.L.) to the data, which for 6 — 1 = 5 degrees of freedom, translates 
to x 2 /dof < 2.2. For each operator, we find upper bounds on A?- as given in Ta- 
ble |l]. OPAL has recently announced 95% C.L. lower bounds on the scales A of each 

1 Thus our remarks do not apply to proposals such as that of Ref. [||. 
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Table 1: 95% C.L. upper bounds on the A?- in TeV. Dashed entries are operators 
which do not fit the HERA data for any A. 



of these operators ||. Given these results, only the A^ R and A RL operators (with 
f] = +1) can have sufficient strength to explain the HERA data, and simultaneously 
satisfy the OPAL bounds. Of course, neither allowed operator will produce a peaked 
distribution in x. 

In any case, there is another strong constraint on these interactions, arising from 
atomic parity violation (APV) measurements. Experiments using cesium (Z = 55, 
N ~ 78) have found its weak charge to be Q^ p (Cs) = —71.04 ± 1.81, while the 
SM predicts Q^(Cs) = -73.12 ± 0.09 JK|. A contact interaction shifts Q w by an 
amount [fL4l: 



AQ w = -2[C lu (2Z + N)+C ld (2N + Z)\ (2) 

where 

r V2tt / 7] q RL _ rf LL _ r] g LR rf RR \ 

lq G F {(A RL r (AL) 2 (AL) 2 (A^)V' [) 

For any one operator with A < 3 TeV, AQw ~ ±20, which is decisively excluded 
by the APV experiments. Logically, we cannot exclude the possibility of a conspir- 
acy, with cancellations among individually large contributions from several operators 
(this can be natural however for contact interactions generated by a parity-conserving 
gauge interaction), but clearly these contact terms provide at best a marginally con- 
sistent parameterization for understanding the HERA data. 

A w-channel leptoquark $, in the heavy mass limit, goes over into the same types 
of contact interactions just discussed. But in this case an additional problem for 
light leptoquarks presents itself within the HERA data. Given e + p scattering in 
the li-channel with a valence (or sea) quark, by crossing and charge-conjugating the 
Feynman diagram, one arrives at e~p scattering in the s-channel off the same valence 
(or sea) quark. Therefore a(e~p) ~ (M$/r$)a(e + p) due to resonant production of $ 
in e~p scattering, where T$ is the leptoquark width: 

A 2 A 2 

Tscalar LQ = n^-M L Q, Tvcctor LQ = ^^j^^ 1 ^' ^ 

where n is the number of available leptoquark decay channels. Given typical values 
T$ ~ 10 MeV, then even with its much smaller integrated luminosity, this represents 
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Operator 



SU(3) SU(2) U(l) 



CC mode 



3 1 1/3 

3 3 1/3 

3 2 7/6 

3 2 1/6 

3 2 -7/6 

3 1 -1/3 

3 1 -4/3 



Yes e 

Yes e~ 

No e+ 

No e+ 

No e+ 

No e~ 

No e~ 



Table 2: List of scalar leptoquark operators. For each operator, the S77(3) x SU(2) x 
C/(l)y quantum numbers of the leptoquark, $, are shown. The fifth column indicates 
whether HERA should find CC events, and the final column lists the mode in which 
HERA should dominantly produce the given leptoquark. Q and L represent 577(2) 
doublet quarks and leptons, while e, u and d are SU (2) singlets. 

a signal of about 10 3 events in e~p in the current data. Thus u-channel leptoquark 
exchange cannot be responsible for the HERA anomaly. 

Thus s-channel production is the most attractive option. It is, moreover, con- 
sistent with the peaked distribution in x observed at HERA. In this channel an 
important question arises, whether the e + p scattering is off valence or sea quarks. If 
the scattering were off sea quarks, then by charge- conjugating the relevant diagram, 
one would arrive at e~p scattering off valence quarks, also resonant. At x ~ 0.5, 
the valence and sea parton distribution functions come in the approximate ratio 
u v « 4d v « 200w s = 200(7,. (Recall that as x —> 1, u v (x) « 2d v (x)/(l — x).) So 
for every one anomalous event observed in the e + p mode, HI and ZEUS should have 
found between 2 to 9 such events in the e~p mode (depending on whether the coupling 
is to u or c?-quarks), despite the small integrated luminosity. Since no such excess 
has been observed, the e + p anomaly must be described by s-channel production of a 
leptoquark coupling to valence quarks in the proton. 

Scalar Leptoquarks. For scalar leptoquarks, there are only a few operators which 
are gauge- invariant, renormalizable and couple to both electrons and quarks. They 
are listed in Table |2|. Also shown are the quantum numbers of the leptoquark, $, 
under S77(3) x SU(2) x U(l)y. In the column "CC" is indicated whether or not 
the operator leads to CC events at HERA in addition to the NC. In the last column 
("mode") is listed the dominant production mode for the given leptoquark. We 
already argued that the events at HERA must be coming from s-channel scattering 
off of valence quarks in the proton; such operators are listed as "e + " mode in the 
table. The remaining "e~" operators would have been already detected in HERA's 
e~p data given the size of coupling necessary to explain the anomaly in the e + p data. 

Note that if a leptoquark couples significantly to both LH and RH quarks, it gen- 
erates effective 4- Fermi operators of the form W^difi^y^. These lead, for example, 
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to helicity-unsuppressed tt — > ev decays, which are severely constrained by exper- 
iment ||14|| . Among the operators relevant at HERA, this means that we may not 
identify <£> e Q = &* Lu (where the subscripts indicate the SM fermions to which the 
leptoquark couples ) despite their identical quantum numbers. Alternatively, for a 
$(3, 2, 7/6) we must forbid one or the other of its Lu$ or eQ§* couplings. 

Of the 7 leptoquark operators listed in Table |2|, only 3 are consistent with the 
HERA data as we interpret it: & Ld and $ 6 q. All three are SU{2) doublets, 

and we will assume throughout the rest of this analysis that the member coupling 
to electrons is the lighter so that it does not decay to the component coupling to 
neutrinos and a real or virtual W . Such a decay would have a very different signature 
at HERA. We also note that the leptoquark which can appear in a 5 of SU (5) is not 
among the three possible states useful for explaining the HERA data; however, the 
$m leptoquark can fit into a 10 of SU(5). 

Notice that all scalar leptoquarks which would have given a CC signal have been 
ruled out already by lack of NC events in e~p mode. Thus for scalar leptoquarks, we 
conclude that there will be no CC signal at HERA. Conversely, if a CC signal is seen 
at HERA, a scalar leptoquark interpretation is strongly disfavored. 

We performed a similar \ 2 analysis here to the one described above for contact 
interactions, except there are now 2 independent variables: Mlq and giQ. Relevant 



cross-section calculations have been performed in Ref. [|ll|]. Here again the APV 
measurement can strongly constrain the range of parameters allowed. For scalar 
leptoquarks one finds 

For the leptoquarks of interest, the values of (5z,8n) are: (-2,-1) for $ e g, (2,1) for 
§Lu, and (1,2) for & Ld . For a given mass and coupling, the constraint on $ e g will be 
weaker than for the others, since it shifts Qw in the direction of experiment. 

Given the 6 + 1 — 2 = 5 degrees of freedom, we explore the range of (M L q, g L q) 
which give \ 2 jdoj < 2.2 (i.e., 95% C.L.). The 95% confidence regions are shown in 
Figure |I]for each of <3> e Q, &l u , and For &l u and §Ld, the leptoquark couplings are 
constrained to be small by the APV bounds discussed above. For $ e g, we assume that 
the members of the leptoquark doublet are nearly degenerate and both are produced 
at HERA; if one or the other is significantly heavier, then the best fit values as far as 
HERA is concerned go over to those of $l„ or $Ld, as appropriate. It may seem that 
in the combined HERA + APV analysis for $ e Q, large values of the corresponding 
coupling, g L Q, still provide good fits since they improve agreement of theory with 
the APV data. However, as we discuss later, a combination of rare decay constraints 
coming from Kl — > it°£ + £~ and D° — > tt°£ + £~ exclude values of g^Q > 0.08 for this 
operator for masses in the interesting range. 

We also calculate "la" bounds on Mlq and giQ by demanding x 2 < Xmin + 1- 
These bounds are not true la bounds since the iso-x 2 lines are far from elliptical but 
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Figure 1: 95% CL. regions in the (Mlq, Qlq) plane derived from fitting to the HERA 
Q 2 distribution and APV in cesium. The three scalar leptoquark operators are shown: 
$l„ (solid), $ e Q (dotted), and &Ld (dashed). For $ iu and $m the best fit values 
for Mlq and qlq are shown with la error bars. For $ e g, couplings above 0.08 are 
excluded by a combination of rare decay constraints. 

we consider the bound indicative of where the best fit values fall. For the &lu and 
§Ld leptoquarks, one finds: 

M LQ = (200 ± 8) GeV, 9lq = 0.025 ± 0.005 ( ] 

®Ld- Mlq — (195l§) GeV, qlq — 0.05 ± 0.005 W 

A similar bound on the $ e Q leptoquark is difficult to determine. The \ 2 is small 
over a wide range of masses (190 GeV < Mlq < 260 GeV) and couplings (0.01 < 
Qlq ^ 0.08), and once again it is only the x-distribution at HERA which can further 
constrain them. 

It is noteworthy that the fits have no information put into them a priori about 
the observed distributions in x. Nevertheless, for the and $^ leptoquarks, the 
best fit values to the shape of the Q 2 distribution are at Mlq consistent with the 
Mlq = \/~xs extracted from the HERA data. We consider this an additional piece of 
evidence in favor of an s-channel interpretation of the HERA results. 
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Vector Leptoquarks. The HERA data, considered by itself, can be similarly well 
described by an s-channel exchange of a vector leptoquark with Mlq ~ 200 GeV 
and coupling smaller than its corresponding scalar leptoquark. However, such 
states are disfavored by an analysis of Run I data taken by the DO collaboration 
at Fermilab, as we now argue. First, note that DO has reported a bound on 1st 
generation scalar leptoquarks of 175 GeV at the 95% confidence level [^]. This 
corresponds to a bound on the leptoquark production cross-section (including the 
effect of cuts and efficiencies) of a < 0.4 pb Now, the production cross-section 
for a pair of vector leptoquarks, $ A1 , at the Tevatron has a model dependence that 
stems from the possibility of a coupling n^G^ u ^ v to the gluon field strength tensor 
G^ u . (There is an insignificant contribution to the production cross section from the 
leptoquark Yukawa couplings given the determination of their size from the HERA 
data.) It can be seen from the analysis of Hewett, et al [JT5 j that the contribution from 
k destructively interferes most strongly with the process qq — > g —>■ $'$ at values 
k ~ —0.5. For this most conservative value, the production cross-section for pairs of 



vector leptoquarks at the Tevatron exceeds 0.4 pb for masses M L q < 215 GeV [[15 
Thus such masses are excluded at the 95% confidence level. If on the other hand 
we take the value of k to lie in the range (0.0-1.0), then masses Mlq < 240 GeV 
are excluded. We therefore view vector leptoquarks as only a marginally consistent 
possibility at best. 



Flavor Violation. The three effective interactions of scalar leptoquarks that we 
have identified as capable of explaining the HERA data are 

(Lowercase Q 's correspond to the components of the leptoquark SU{2) doublets of 
electric charge Q.) Here the operators are written in the mass eigenbasis, where i,j 
are the generation indices and V in O" is the CKM matrix. The flavor structure of 
these Yukawa matrices (A, A', A") is severely constrained by limits from rare processes 
(for earlier work in this regard, see ]T4{|). 

To be definite take the operator O first. If O explains the HERA data, An ~ 0.05 
is needed. The other elements of A are then constrained as follows. The decay 
Kl — > e + e~ occurs at tree-level by the exchange of the (j) 2 / 3 leptoquark. However, 
the rate is strongly suppressed by a helicity factor (rn 2 /m 2 K ) and therefore does not 
give any useful limit. On the other hand, the decay Ki — > n°e + e~, which arises 
from the same tree-level quark diagram, is helicity unsuppressed. The limit on this 
decay, BriK^ — > 7r°e + e~) < 4.3 x 10~ 9 , translates into a bound of A11A12 < 1.5 x 
lO- 4 (M 02/3 /2OOGeV) 2 . For An w 0.05, this means \ 12 < 3 x lO" 3 . Such a hierarchy 
pattern seems quite unusual; it deviates considerably from the pattern of Yukawa 
couplings that we have become used to in the quark sector. 



7 



The doublet partner of 2//3 must be nearly degenerate with it, else their 

mass-splitting will contribute to the p-parameter excessively. The tree-level exchange 
of causes the decay K + — ► ix + vV. The experimental limit Br(K + — > -n + uV) < 

2.4 x 10" 9 leads to the bound A;iA 2i < 1.2 x Kr 4 (M^-i/ 3 /200 GeV) 2 . One again sees 
a severe limit on the second generation couplings. 

The leptoquarks (p 2 ^ 3 and mediate the process \i — > e + 7 through one loop 

induced magnetic moment interactions. The present limit on the branching ratio, 
Br(ji -> e + 7) < 4.9 x KT 11 , implies that A u A 2 i < 1.6 x 10~ 4 (M$/200 GeV) 2 . 
Similarly, the decay Kl — > /ie can occur at tree level if both An and A22 are nonzero. 
The limit on this decay, Br(Ki — > /xe < 3.3 x 1CT 11 ), translates into a bound of 
A11A22 < 2.8 x 10~ 6 (M$/200 GeV) 2 . Evidently it is not enough that the couplings are 
flavor diagonal; individual quark and lepton number must be preserved. 

There are analogous constraints on the coupling AL of O' coming from D° — > 
7r°e + e~, etc, but they are weaker numerically, by about a factor of 5, in the Yukawa 
couplings. The couplings X"j in O" have the interesting feature that one cannot 
simultaneously suppress D° — > 7r°e + e~ and Kl — > 7T°e + e~. Suppressing the more 
strongly constrained mode Kl — > 7r°e + e~ leads to a limit on the individual diagonal 
Yukawa coupling |A^| 2 < 4 x 10 _3 (M$»/200 GeV) 2 from D° -> 7r e + e-. This is the 
bound imposed on the HERA + APV fit in Figure [I] for the $ e g operator. There 
is also a stringent limit on the product (X" T V)n(X" T V)l 2 from Kl — > n°e + e~ decay. 
In particular, if only A'^ is nonzero (in this basis), the constraint \X"i\ 2 < 6.8 x 
10~ 4 (M$///200 GeV) 2 , will nearly exclude the possibility of explaining the HERA 
data with this operator. 

Additional constraints arise from box diagrams which induce K° — K mixing, 
D° — D mixing, etc, but these are less severe than the ones quoted above, essentially 
because they scale as A 4 in the amplitude. One process worth noting is the box 
diagram contribution to Kl — > Tx°i + C.~ and D° — > tt £ + £~ with the exchange of one 
W and one leptoquark. Such a diagram directly constrains the diagonal leptoquark 
coupling, although the limits are rather weak: | A 2 2 1 2 < 0.06(M$/200 GeV) 2 from 
Kl -> 



Conclusion: A Model. Limits arising from flavor physics provide very severe con- 
straints on any model containing leptoquark degrees of freedom capable of explaining 
the HERA data. In the absence of a general principle, it strains credulity to pos- 
tulate the required large number of tiny parameters in a fundamental theory. An 
adequate interpretation of the HERA data must address this question. Neither the 
direct introduction of fundamental leptoquark degrees of freedom, nor other possibil- 
ities such as i?-parity violating couplings of MSSM squarks [0, naturally explains 



this remarkable flavor structure; namely that the leptoquarks must, to a very high 
degree of accuracy, couple diagonally to the mass^eigenstate quarks and leptons, and 
conserve individual quark and lepton numbers. Clearly what is required is that the 
underlying physics responsible for the HERA anomaly is universal in character with 
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regard to the generations. 

The natural implementation of this universality is a gauge principle, with the 
three generations having identical charge assignments. We have seen that the HERA 
phenomenology disfavors t-channel Z' exchange as an explanation compared to the 
production of an s-channel resonance. We are therefore led to consider the possibility 
that a new strong short-ranged gauge interaction leads to the formation of relatively 
light leptoquark quasi-bound states. Moreover, if the unification of couplings observed 
in the MSSM is to have any chance of being maintained, this gauge interaction must 
commute with the Standard Model interactions. As an example consider an additional 
U(l) interaction at a strong coupling fixed point, spontaneously broken at the TeV 
scale. While one cannot reliably calculate the consequences of such an interaction, 
a not unreasonable hypothesis for the dynamics of such a theory might be that the 
low energy spectrum is saturated by resonances having small residual interactions. If 
the original gauge interaction commutes with 577(5), then the resonances automati- 
cally fill out several complete SU(5) multiplets, each having different flavor quantum 
numbers. More importantly, the generation independence of the U(l) gauge coupling 
implies that the masses and couplings of the resonances (to free quarks and leptons) 
are to a good approximation invariant under an SU(3) symmetry that rotates the 
generations. This symmetry is only broken by effects that depend on the usual Higgs 
Yukawa couplings, h, of the SM as h 2 . The flavor physics limits discussed above are 
easily and naturally satisfied. 

A mo del- independent prediction emerging from this scenario is the existence of 
a large multiplicity of very nearly degenerate resonances. This is to be expected 
from any model which preserves so carefully the 577(3) flavor symmetries among the 
generations, as seems to be required by the experimental flavor constraints. These 
states could be discovered, or excluded, in future experiments at the Tevatron or 
LEP II. There is already an interesting constraint arising from FNAL data, since 
at minimum there are two degenerate states containing an electron in combination 
with either a 1st or 2nd generation quark. From the DO bound on leptoquarks which 



couple to electrons |T2[ the total production cross section to these states has to be 
less than 0.4 pb, which translates to a limit of M^q > 190 GeV. Degenerate with 
these states will be leptoquarks composed of \iq and rq, although the present bounds 
on such states are much weaker. 

Note added. While this paper was in preparation, several other papers inspired by 
the high-Q 2 HERA anomaly appeared [17]. Compared to other authors, we appear 



to place greater emphasis on the unusual flavor structure required. 
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